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Abstract 

In this paper, we consider the problem of maximizing the energy efficiency (EE) for multi-input 
multi-output (MIMO) interference channels, subject to the per-link power constraint. To avoid extensive 
information exchange among all links, the optimization problem is formulated as a noncooperative game, 
where each link maximizes its own EE. We show that this game always admits a Nash equilibrium 
(NE) and the sufficient condition for the uniqueness of the NE is derived for the case of large enough 
maximum transmit power constraint. To reach the NE of this game, we develop a totally distributed EE 
algorithm, in which each link updates its own transmit covariance matrix in a completely distributed and 
asynchronous way: Some players may update their solutions more frequently than others or even use 
the outdated interference information. The sufficient conditions that guarantee the global convergence of 
the proposed algorithm to the NE of the game have been given as well. We also study the impact of the 
circuit power consumption on the sum-EE performance of the proposed algorithm in the case when the 
links are separated sufficiently far away. Moreover, the tradeoff between the sum-EE and the sum-spectral 
efficiency (SE) is investigated with the proposed algorithm under two special cases: 1) low transmit 
power constraint regime; 2) high transmit power constraint regime. Einally, extensive simulations are 
conducted to evaluate the impact of various system parameters on the system performance. 
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1. Introduction 

Past few years have witnessed tremendous advaneement in wireless eommunieations, ineluding 
the significant improvement of transmission rate [[T|, Q. However, the impact of power con¬ 
sumption on the environment is neglected. It is reported that the total energy consumption of the 
communications takes up more than 3 percent of the worldwide electric energy consumption [[^ 
and the portion is expected to increase due to the explosive growth of high-data-rate applications 
in the future. Hence, energy efficiency (EE) has gained lots of attention and will be one of key 
issues in future fifth-generation (5G) mobile networks Q. On the other hand, the interference 
channel (IC) has been modeled mathematically for many practical systems where multiple 
uncoordinated links share the same channel, such as femtocells, ad hoc wireless networks, 
cognitive radio, etc Q, Q. Eurthermore, due to the development of advanced multi-antenna 
techniques Q, each transmission node is able to accommodate multiple antennas 0, 0. It 
is well known that multi-input multi-output (MIMO) system has the great potential for providing 
high SE by employing spatial multiplexing techniques p0| . Hence, it is of great importance to 
study the energy efficient transmission strategy in MIMO ICs. 

This paper focuses on the EE maximization problem for MIMO ICs with per-link transmit 
power constraint. To solve this problem, one may consider centralized solutions, which require a 
central processing unit (CPU) to collect all complex-valued channel matrices over the network. 
The CPU will compute all links’ transmit covariance matrices and send them to the corresponding 
links. Hence, for large-scale networks, the centralized approaches suffer from heavy feedback 
overhead and high computational complexity, which hinders practical implementations. More¬ 
over, there may not exist a CPU for some wireless networks, such as ad hoc or wireless sensor 
networks. 


Recently, distributed algorithms to deal with this problem attract intensive attentions [ 111-[ 14|. 
Here, “distributed” means that precoders can be computed at the transmitters with only local 
channel knowledge and limited (or no) information exchange over different links. Generally, 
distributed processing for MIMO systems has the benefits of low communication exchange 
overhead, low computational complexity, more scalability, low system costs, etc. The classical 
distributed algorithm based on dual decomposition technique is designed to decompose the 
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coupling constraints among links O, [[Tg. Specifically, by introducing the dual variables 
assoeiated with the eoupling eonstraints, the original problem ean be divided into several in¬ 
dependent subproblems, eaeh of whieh ean be solved in a distributed way. Then, all the links 
exehange some neeessary information to update the dual variables. For more praetieal networks 


with individual link power budget eonstraints, [15| devised a deeentralized beamforming EE 
(DBEEE) algorithm for symmetrie MIMO ICs, where the distanee from one transmitter to its 
desired reeeiver is identieal for all links, and eaeh transmitter has the same distanee to all its 
unintended reeeivers. In eaeh iteration, all the reeeivers should feed baek the equivalent ehannel 


matriees to all the transmitters in the network. In [I6|, the authors designed a two-layer EE 
(TEEE) algorithm based on the generalized weighted minimum mean square error (WMMSE) 


approaeh p4| : The inner layer to update preeoders/deeoders; the outer layer to update some 
parameters. Similar to the algorithm in p3| , in eaeh inner iteration eaeh reeeiver should feed 
baek the updated weight matrix and the positive definite eovarianee matrix to all the transmitters 
in the network. However, for these distributed algorithms, in eaeh iteration eaeh reeeiver needs 
to eompute the neeessary eomplex-valued matrix and then feed it baek to all the transmitters in 
the network, whieh eould induee serious implementation ehallenges sueh as a large amount of 
feedbaek overhead, poor sealability and heavy eomputational burden at the reeeivers. Moreover, 
all the links should be synehronous, whieh is diffieult to be satisfied, espeeially for large-seale 
ad hoe networks or wireless sensor networks. One novel distributed algorithm based on the 
adaptive priee was proposed in to deal with the weighted sum EE maximization problem 
for single-input single-output (SISO) ICs. 

Henee, one distributed algorithm with mueh lower feedbaek overhead is more desirable. 
Noneooperative game theoretieal approaeh has been reeognized as a powerful tool to devise 
totally distributed algorithms, in whieh eaeh link just maximizes its utility without the need of 
information exehanges among the links. A number of researehes have applied game theory to 
design energy effieient eommunieations for ICs p^-[|2T| or multiple aeeess ehannel (MAC) 

The EE optimization problem in flat fading single-input single-output (SISO) ICs was eonsidered 


in [ 191, where one distributed algorithm based on non-eooperative game was proposed. Both the 
existenee and uniqueness of Nash equilibrium (NE) were analyzed. This work was extended to 
frequeney-seleetive ehannel in p0| and to a relay ehannel in [j^. Although p0| proved the 


existenee of the NE, regarding the uniqueness of the NE, [20| only showed that the number of 
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NEs is determined by the eross-channel gains and the direet ehannel gains, without quantifying 
how they are related with eaeh other. All these studies eonsidered single-antenna ICs and they 
apply for multi-antenna ease only if the transmit powers are optimized with fixed transmit 


direetions as shown in |231. In [22|, the authors considered the EE maximization problem for the 


MIMO uplink systems with each user transmitting only one stream. This problem is formulated 
as a non-cooperative game, where the uniqueness of NE is guaranteed by the fact that the EE 
function is S-shaped. 

In this paper, we apply the non-cooperative game theoretical approach to deal with the EE 
maximization problem for MIMO ICs, where each link attempts to maximize its own EE by 
jointly optimizing transmit power and beamformers. It is a nontrivial extension of the SISO case 


in [20|. There is an explicit relationship between the power allocation among different links 
and the achievable rates in the SISO case. This property is critical in deriving the conditions 
of the existence and uniqueness of the NE by using the standard function p0| . However, in 
MIMO systems, this relationship is implicit, as power allocation is carried out through matrix 
manipulations. Moreover especially when multiplexing is utilized with the MIMO, different from 


the beamforming case in [22| for MAC, the EE design is shown more general for the MIMO IC 
and hence more difficult since there exist both inter-node (mutual) interference and intra-node 
inference in the MIMO IC using multiplexing. 


A. Related Work 

Recently, distributed algorithms for MIMO ICs have been extensively studied in the literature, 
such as linear iterative approximation (EIA) algorithm pA| , the WMMSE algorithm [[M|, nonco¬ 
operative game theoretic algorithm p^ , etc. Eor the EIA algorithm, it is designed based on the 
first-order Taylor expansion of the non-convex part of the weighted sum spectrum efficiency (SE) 
objective function. Sequential convex optimization approaches were then presented for dealing 
with various scenarios, e.g., the multi-band scenario in p5|, the MAC in [[26|, the broadcast 


channel in [271, the cognitive radio networks in [ 12|. However, in the EIA algorithm only one user 
is allowed to update its covariance matrix at one time, which may lead to significant latency 
especially in dense networks. By establishing the equivalence between the weighted sum SE 
problem and weighted sum mean square error minimization problem, [ [T4| proposed the WMMSE 
algorithm that allows multiple users to update simultaneously. In this algorithm, the local optimal 
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solution is obtained via alternatively optimizing the linear transceivers and iteratively updating 
the weight matrices. The authors show that when the utility function satisfies some conditions, 
the algorithm is guaranteed to converge to the stationary point of the original problem. The 
WMMSE algorithm has been applied in various setups, please see [ |28| and references therein. 

However, to successfully implement the ILA and WMMSE algorithms in a distributed manner, 
two assumptions are required: 1) Perfect channel reciprocal between the uplink phase and the 
downlink phase (in time-division duplexing mode); 2) Synchronization between all the links. 
In practice, the communication systems usually operate in frequency-division duplexing (EDD) 
mode. The channel reciprocal is thus hard to achieve. In addition, for wireless senor or ad hoc 
networks, asynchronous among the links is more desirable. More importantly, in each iteration 
each receiver needs to calculate the pricing matrix in the lEA algorithm or the weight matrix in 
the WMMSE algorithm, and then feeds them back to all the transmitters in the network, which 
posses serious implementation issues such as a large amount of heavy feedback overhead, poor 
scalability and heavy computational burden at the receivers. 

On the other hand, non-cooperative game theoretical approaches have attracted extensive 
attentions. Studies in this direction are plentiful in literature, e.g., [ [29| for the MIMO MAC, 
[2^, [[^-[[^ for the MIMO ICs. In p9|, the authors formulated the problem as a non¬ 


cooperative game and the authors proved that each user’s optimal eigenvectors do not depend on 
the channels of others. Based on this fact, the authors showed that the existence and uniqueness 
of NE is guaranteed when the numbers of transmit and receive antennas become large. The 


main technique in [29| is random matrix theory. Eor the MIMO ICs, the authors in [24| first 
utilized the noncooperative game framework to deal with the SE maximization problem, where 
the iterative water-filling algorithm was proposed to find the NE of the game. However, the 
existence of the NE was only shown by the simulation results without theoretic guarantees. 


neither the NE uniqueness. Then, in [30|, Scutati et al proved the existence of NE and provided 
sufficient conditions for the uniqueness of the NE, which can be checked in practice. However, 
the results are only valid for square nonsingular channel matrices. Eater on, they generalized the 


results to a more general case with arbitrary channel dimensions in [311. The cognitive radio 


network with null shaping constraints on the primary user is considered in [32| and its robust 


version in [331. Most recently, these works were extended to the multicell case in [34|, 


with multiple users per cell. In [36 j, the authors formulated the SE maximization problem as a 
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cooperative game. Speeifieally, by fixing the outgoing eooperative set and ineoming eooperative 
set, the authors formulated this problem as a non-eooperative game, where the existenee and 
uniqueness of the NE was analyzed. Then, the eoalitional game theory [|T7| was applied to obtain 


the stable of the eooperative set. The work that is most elosely related to ours is [311, where the 
noneooperative game was formulated for the SE maximization problem of the MIMO ICs and 
one asynehronous distributed algorithm was proposed to reaeh the NE of the game. 

In eontrast to the most of the above eited papers whieh foeus on the (weighted) sum SE 
problem, in this work we eonsider the EE maximization problem. Eor SE optimization problems, 
it is known that all the transmitters use full power during transmission in order to maximize its 
own SE. Based on this faet, the best response strategy at the NE ean be written in a elosed-form 
water-filling solution, whieh ean be interpreted as a projeetor on the eonvex and elosed set. This 


interpretation enables the authors to derive the uniqueness of the game’s NE [31 j. However, the 
study of EE maximization problem eannot be obtained by employing the methodologies sinee 
the transmitters in faet use a portion of the power, instead of full, to aehieve energy effieient 
transmission. 


B. Contributions 

In this paper, we apply the non-eooperative game theoretieal approaeh to deal with the EE 
maximization problem for MIMO ICs, where eaeh link attempts to maximize its own EE by 
jointly optimizing transmit power and beamformers. 

The main eontributions and observations of our work are summarized as follows. 

1) The EE maximization problem in MIMO interferenee ehannels is modeled as a noneoop¬ 
erative game where eaeh MIMO link eompetes against the others by ehoosing its transmit 
eovarianee matrix to maximize its own EE. We show that the NE of this game always 
exists and derive suffieient eonditions for the uniqueness of the NE for the ease of large 
enough maximum transmit power eonstraint. 

2) To reaeh the NE of the game, we provide a totally distributed EE algorithm named 
Asynehronous Distributed Energy-Effieient (ADEE) algorithm, whieh is the extended ver¬ 
sion of simultaneous updating proposed in [ [20| . In this algorithm, all users apply the 
fraetional programming to update the transmit eovarianee matriees and these updates ean 
be performed in a totally asynehronous way, whieh means some links may update their 
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transmit covariance matrices more frequently than the others and they may even use the 
outdated information of the measurement of the interferenee generated by the other links. 
In addition, during the updating proeedure of the algorithm, there is no need for the links to 
exehange the signaling overhead mutually. These features make our distributed algorithm 
more appealing for praetieal implementations. We provide the suffieient eonditions for the 
global eonvergenee of this algorithm to the unique NE of the game. Interestingly, we find 
that these eonditions eoineide with the eonditions for the uniqueness of NE. 

3) We study the impaet of the eireuit power eonsumption on the overall SE and EE perfor- 
manee of the system for one speeial ease when the links are separated suffieiently far away. 
We show that the overall SE inereases with the eireuit power eonsumption, but the overall 
EE deereases with it. Although this trend is derived for this speeial ease, from simulations 
we find the trend holds for the general ease when the interferenee among the links is 
sizeable. This observation implies that when the eireuit power eonsumption inereases, we 
should enhanee the transmit rate or SE in order to obtain the best EE performanee. 

4) The tradeoff between SE and EE is investigated for the proposed algorithm (denoted as 
ADEE algorithm) and the SE maximization algorithm (denoted as ADSE algorithm) in 
10. Two speeial eases are studied: the transmit power eonstraint approaehes zero or 
infinity. Eor the ease of low transmit power eonstraint, we show that both algorithms use 
all power to transmit and thus aehieve the same performanee in terms of the overall SE 
and EE performanee. However, for the latter ease, the ADSE algorithm always uses all 
available power to transmit, yielding severe inferenee over the network. In this ease, the SE 
aehieved by the ADSE algorithm will not inerease. Then, the EE aehieved by the ADSE 
algorithm will approaeh zero due to the signifieant power eonsumption. On the other hand, 
for the EE metrie, the ADEE algorithm is unwilling to eonsume all power in this ease. 
As a result, the SE and EE aehieved by the ADEE algorithm will beeome eonstant in the 
ease of the high transmit power eonstraint. 

The rest of the paper is organized as follows. In Seetion |I^ we introduee the system model 
and formulate the optimization problem as a strategie noneooperative game. Then, we show 
that this game always admits a NE and derive suffieient eonditions for the uniqueness of the 


NE in Seetion III To reaeh the NE, a totally asynehronous and distributed algorithm is given 



in Section IV In Section |V| we study the impact of the circuit power consumption on the 
system performance in terms of the sum-SE and sum-EE, along with the study for the tradeoff 
between the sum-EE and sum-SE for the proposed algorithm. Section |V^ provides representative 
numerical results to study the effects of different system parameters on the proposed algorithm. 


Einally, some conclusions are drawn in Section VII 

Notations: (•)*, (•)^, (•)^, vec(-), E{-}, tr() and ® are conjugate, transpose, Hermitian 
transpose, stacking vectorization operator, expectation operator, trace operators and the Kronecker 
product operator [ |38| , respectively. Uppercase and lowercase boldface denote matrices and 
vectors, respectively. Eor matrix A, [A],^ and [A]jj represent the column of A and the 
(i,j) element of matrix A, respectively. ||A ||2 denotes the spectral norm of A [38|. ||A||p, 
denotes the Erobenius norm of A. Amin(A) stands for the minimum eigenvalue of A. A ^ B 


means A —B is positive semidefinite. The spectral radius of A is denoted by p(A) [381. rank(A) 
denotes the rank of A. Eor vector a G ||a ||2 represents the Euclidean norm defined as 

||a ||2 = Va-^a. The sets of m x n complex matrices, n x n complex positive semidefinite and 


definite matrices are denoted by and respectively. DxY denotes the Jacobian 


matrix of function Y with respect to (w.r.t.) X [39|. I and 0 represent the identity and zero 
matrices with appropriate dimensions, respectively, [x]’'" is equivalent to max{0,x}. 


II. System Model and Problem Eormulation 

Consider a A-link MIMO interference channel with K transmitter-receiver pairs. Each link 
consists of one transmitter with M transmit antennas and one receiver with N receiver antennas. 
All links are simultaneously communicating over the same channel. At receiver k the received 
complex baseband signal vector G is given by 

Eh jkXj -h rifc, (1) 

where x^. G denotes the transmit signal vector of link k, Ylkk ^ is the direct channel 

matrix of link k, G is the cross-channel matrix from transmitter j to receiver k, and 

rifc G is circularly symmetric, zero-mean, complex Gaussian noise with normalized identity 


'“Receiver k” represents the receiver of the fc"* link. In the following, “transmitter fc” means the transmitter of the link. 
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covariance matrix. For each link, the total average transmit power should satisfy the per-link 
power eonstraint: 

Pk = tr(Qfc) < Pt, (2) 

where Qk = E {^k^k} is the covariance matrix of and Pt is the maximum transmit power. 

To reduee the eomplexity of deeoding at the reeeivers, it is assumed that joint deeoding of 
the interfering signals is not an option and the interferenee is treated as noise at the reeeivers. 
Thus, the SE of link k is given by (in bit/s/Hz) 


Cfe(Qfc, Q-fc) — log2 |I + ^Hk^kQkl, 


(3) 


where Rfc = l + represents the interferenee-plus-noise (IPN) eovarianee matrix 

at reeeiver k and Q_fc = [Qi, • • ■ , Q^-i, Qfc+i, ■ ■ ■ , Qk] denotes the set of all links’ eovarianee 
matriees, exeept that of link k. In this work, it is assumed that eaeh reeeiver k ean perfeetly 
measure the IPN eovarianee matrix R^ and estimate the direet ehannel and then report 
them baek to transmitter k. The ehannels are assumed to vary suffieiently slowly sueh that it 
ean be eonsidered fixed during the transmission. 

In order to design energy effieient transmissions, the total power eonsumption should be 
eonsidered at eaeh link k: power used for reliable data transmission P*., eireuit power during 
transmission Pc, whieh is the power eonsumed by the mixers, filters and digital-to-analog 
eonverters, digital signal proeessing (DSP), ete. The power eonsumption of the DSP depends 
on the number of eomputations of the algorithm and the signaling overhead. It is diffieult to 
aeeurately model this kind of power eonsumption. Thus, the eireuit power Pc is modeled as a 


eonstant in this work for the sake of analysis, as simplified in most of the existing works [15|, 


|19|-[21|, [231. Even though, we will evaluate the system performanee under different Pc via 
simulations. Then, the EE (in bits/Hz/Joule) of link k, defined as the ratio of SE to the total 
power eonsumption, is given by 

CkiQk, Q-fc) 


EEfc(Qfc, Q-fc) — 


(4) 


Pfc+Pc 

Sinee our goal is to devise a totally distributed algorithm that requires neither a CPU nor 
information exehange among the links, we formulate the optimization problem as the following 
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noncooperative game: 


max EEfc(Qfc,Q_fc), 

(Q) : Qfc Vfc G ^K, 

s.t. Qfc G Wfc 


( 5 ) 


where EEfc(QA;, Q-fc) is the payoff funetion of link k defined in (4), ^k = {1, 2, • • ■ , K} is the 
set of all links and Wk is the set of admissible strategies of link k, defined as 


m = {Qk e e "'^,tr(Qfc) < Pt}. 


( 6 ) 


In this game, eaeh user eompetes against the others by ehoosing his own eovarianee matrix 
that maximizes his own payoff funetion subjeet to the strategy set. A solution of the game to 
reaeh a NE is when eaeh link, given the strategie profiles of the others, does not get any inerease 
in its objeetive by unilaterally ehanging its own strategy and is formally defined as follows. 

Definition 1: A strategie profile Q* = (Qfc)A:e'i'x ^ x ■ ■ ■ x y^K is a NE of game Q if 


EEfc(Q^,Q*fc)>EEfc(Qfc,Q%),VQfc G G ^k- 


(7) 


In the fortheoming seetions, we first show that game Q always admits at least one NE. In 
general, game Q may admit multiple NEs, depending on the level of the interferenee from 


the other links [20|. Then, we study the uniqueness eondition of the NE and provide a totally 
distributed algorithm to reaeh sueh a NE. 


III. Existence and Uniqueness of the NE 

A. Existence of NE 

Whether NE exists depends largely on the properties of the payoff funetion. In the sequel, we 
first study the property of the EE funetion and then eheek the existenee of NE. 

Lemma 1: Given the other links’ strategy Q_fc, the EE funetion of link k, i.e., EEfc(Qfc, Q_fc), 
is quasieoneave in Q^. Eurthermore, if the ehannel matrix ^ is full eolumn rank, EEfc(Qfc, Q_fc) 
is strietly quasieoneave in 

Proof: Please see Appendix □ 

Based on Eemma 1, the existenee of NE is given in the following theorem. 

Theorem 1: Game Q always admits at least one NE for any set of ehannels and transmit 
power eonstraints. 

Proof: Please see Appendix □ 
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B. Uniqueness of NE 

In this subsection, we first study the uniqueness condition for the case of full column-rank 
channel matrices. Then, we extend the results to the more general case without making any 
restrictive assumptions on the channel structure. 

1) Case of full column-rank channel matrices 


For the SE maximization game in [311, all the transmitters use full power during transmission 
in order to maximize its own SE. Based on this fact, the best response strategy at the NE can 
be written in a closed-form water-filling solution, which can be interpreted as a projector on the 
convex and closed set. This interpretation enables the authors to derive the uniqueness condition 
of the game’s NE. However, the study of our EE maximization problem cannot be obtained by 


employing the methodologies developed in [311 since the transmitters do not always, and in fact 


not in most cases, use the full power to transmit in our EE case. Here, we consider the case 
when the maximum transmit power is large enough so that each transmitter only uses portion 
of the maximum power at the NE. 

Before providing the uniqueness condition for the EE maximization game, we first introduce 
some useful intermediate results. Given the other links’ strategy Q-k, the best response of the 
transmit covariance matrix of link k is denoted as 


Ffc(Q_fc) = arg max EEfc(Qfc,Q-fc) 

Qfe6Wj; 


( 8 ) 


where function Ffc(Q_fc): W_fc is a complex matrix-valued function with W-k = 

Wi X ■ ■ ■ X Wfc-i X Wfc+i X • • ■ X Wk- We introduce the following mapping function. 


F(Q) = Fi(Q_i) X • ■ ■ X F;,(Q_;,) : W ^ W 


(9) 


where W = Wi x ■ ■ ■ x Wk and Q = [Qi, ■ ■ ■ , Qi^] G W. Using (9) and Definition 1, the NE 


of the game can now be characterized by the following fixed-point equation: 

F(Q*) = QT 


( 10 ) 


According to Proposition 1.1 of [40|, the contraction mappings have a unique fixed point. It 
is formally stated as the following lemma. 

Lemma 2: Game Q has a unique NE if the mapping F(-) defined in (|^ satisfies: 

||F(Q(^)) - F(Q(2))|| < 11q9) _ q(2)|| 


IF’ 


( 11 ) 
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for any two different and in the elosed set W. □ 

However, direetly using this lemma is diffieult. Instead, we give a suffieient eondition for ( [T^ 
to be satisfied, whieh plays a key role in the study of uniqueness of the NE. 


Lemma 3: The suffieient eondition for (12) to hold is that, for any link k, the mapping defined 
in dSl) satisfies 


Fi(QL‘ 1) - Fj(Q 


(2)^ 


< 


K-1 


Q'i - Q 


( 2 ) 




( 12 ) 


Proof: Our main task is to show that, under eondition ( [T^ , the eondition of the mapping F(- 
to be a eontraetion in ([T^ is satisfied. We have 


K 


|F(Q<‘>) - F(Q<^l)||;, = y] ||Fi(Q«) - Ft(QS) 


d2)' 


k=l 


K 


< 




K - 1 


k=l 


= ||q( 1) _ o(2) 




(13) 


(14) 


(15) 


□ 


where ( [T4] ) follows from (12), and ( [T5| ) follows by the definitions of Q_fc and Q. 

Based on the above results, we give a suffieient eondition to guarantee the uniqueness of the 
equilibrium for the ease of full eolumn-rank matriees in the following theorem: 

A / ^ 

Theorem 2: Define T^, = I I + Pt X] '^i,k^fk ) . Suppose that the ehannel matriees 

are full eolumn rank, i.e., rank(Hfc fc) = M. Let ak be 


fcS'I'x 




ak = 


k,k 




, 2 ’ 


(Amin 

where Dq ,,Rfc is the Jaeobian matrix of Rfc w.r.t. Q_fc and is given by 


(16) 


DQ_j.Rfc — ® • ■ ■ , (8) 8) ■ ■ ■ , 8 (17) 

Then for suffieient large maximum transmit power Pt, the NE of game Q is unique if 

ttfc < 


K-1 


,Wk e 


(18) 

□ 


Proof: Please see Appendix 
The importanee of Condition ( [T8| ) is that for given direet ehannel gains, it explieitly qualifies 
how small the multi-link interferenee eaeh link ean tolerate to guarantee the uniqueness of the 


NE. Henee, Condition (18) ean be eheeked in praetiee to faeilitate the admission eontrol. In 
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contrast, [ |20| only showed that the number of the equilibria is determined by the cross-channel 
gains and the direct channel gains, without quantifying how they are related with each other. 

2) Case of more general channel matrices 

In practical systems, the channel matrices may not be full column rank. In this part, we consider 
the more general case without making any restrictive assumptions on the channel structure. 

For each link k, we write the eigendecomposition of 'Hk,k as Hfc^fc=UfcAkV^, where tj^ G 
G are semi-unitary matrices, Ak G is a diagonal matrix with positive 

eigenvalues, and fk is the rank of matrix Ylk,k, i-C-, = rank(Hfc fc) < min{M, A}. To 

maximize each link’s EE in Q, each link fc’s optimal covariance matrix should lie in the subspace 
orthogonal to the null-space of ^ for a given Q_fc. It follows that the best response of the 
transmit covariance matrix of each link k belongs to the following class of matrices: 

Qfc = VfcQfcVf,VA:G^i,, (19) 


with 


e Wfc = {X G G §:^^"^tr(X) < Pt} • 


( 20 ) 


By inserting (19) into game ^ in ([^ and defining fc, game Q can be 

transformed into the following lower-dimensional game Q, given as 


rnax EE^ (Q^, Q-fc) = 
{Q), Q. _ ' _ 

s.t. Qfc G Wfc, 


log2|l+H"^R-^Hfc,fcQfc| 

tr(Qfe)+Pc 


V/c G ^ 


K, 


( 21 ) 


where Q_fc = QiX- • - xQfc.ixQfc+iX- ■ - xQx and = R^^(Q_fc) = 

Now the channel matrices ^k,k £ are full column rank. Hence, we can apply the 

same derivations for Theorem 2 to obtain the uniqueness condition for this more general channel 


case, which is the same as that in (18) except that the channel matrices k are replaced 

by Hj,fcVj, k. 

To give more insights into the physical interpretation of the uniqueness conditions of the 
NE, Fig. is plotted to show the probability for the uniqueness conditions to be satisfied for 
different cases. Specifically, full-rank square, fat and tall channel matrices are simulated. If 


the channel matrices are square or tall, condition (18) in Theorem 2 is applied to check the 


uniqueness condition of the NE. In the case of fat channel matrices, the above derived condition 
for general channel matrices is used to check the uniqueness condition of the NE. The probability 
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of the uniqueness of the NE is defined as the ratio of the number of ehannel matrices that 
guarantee the uniqueness condition to the total number of channel matrices. For simplicity, we 
consider a symmetric system with two links {K = 2) where the direct-channel distances for 
both links are set to be /indirect = 1 and both links have the same cross-channel distances (i.e., 
Di ^2 = D 2 ,i = E*cross)- The power constraint and the noise power are set to be Pt = 10“^W, 
= 10“^W, respectively. Note that Pt ^ a^. These results are obtained by testing over 10000 
random channel matrices whose entries are generated as the circularly symmetric, zero-mean, 
complex Gaussian random variables with variance equal to the square root of the channel path 
loss power, where the path loss exponent is assumed to be 3.5. Several interesting observations 
can be found from Fig.[^ 1) The uniqueness probability of the NE for all cases increase with the 
cross-channel distance corresponding the decrease in the interference. This is reasonable since in 
the extreme case when the cross-channel distance approaches infinity, the system becomes two 
independent point-to-point links. On the other hand, when one of the receivers is too close to one 
of its unintended transmitters, one of the two links should be shut down, which is instructive for 
practical use; 2) For the case of square channel matrices, the uniqueness probability decreases 
with the number of antennas; 3) For the case of fat channel matrices, increasing the number of 
transmitter antennas while keeping the number of receiver antennas fixed leads to an increase 
in the uniqueness probability of the NE. For example, the curve associated to the case of 2 x 5 
MIMO channels is higher than that associated to the case of 2 x 3 MIMO channels for any 
given Across- Similar observations hold for the case of tall channel matrices. This observation is 
of significant importance: equipping unequal number of antennas at transmitters and receivers 
can dramatically improve the uniqueness condition of the NE. 


IV. Asynchronous Distributed Energy Efficient Algorithm 


To reach the NE of game Q, we employ the totally asynchronous algorithm [40|. The main 
characteristic of the asynchronous algorithm is that some users are allowed to update their best 
response more frequently than others. This algorithm has been successfully employed to deal 
with the rate maximization game in 0, 

For the sake of readability, we briefly introduce the asynchronous algorithm and adapt it to 
our EE maximization problem. To this end, we first introduce some definitions and assume that 
the set of times at which the links update their solutions is the discrete set T = {0,1,2, ■ ■ ■ }. 
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cross 


Fig. 1. Probability of the uniqueness of the NE for different numbers of antennas. 


Let C 7~ be the subset of times at which transmitter k updates its solution and let Qfc(f) 
be the updated transmit covariance matrix of transmitter k at time t Denote r^(f) as the 

most recent time at which the interference from transmitter r is measured at receiver k at time 
t. Hence, at time t, link k updates its transmit covariance matrix based on the interference from 


To guarantee that the system is totally asynchronous, three conditions should be satisfied by 
the schedules {T^{t)} and {T^} [40|: Al) 0 < < f; A2) limfc_).ooT,?(4) = +oo; and A3) 

|T^| = oo; where {tk} is a sequence of T^. Based on the above definitions, the asynchronous 
distributed EE algorithm is described in Algorithm 1, where T^ax is the maximum number of 
the iterations. 

Remark 1 - Two special cases: Note that some well-known algorithms such as sequential and 


simultaneous algorithms [40|, where the links update their solutions sequentially and simulta¬ 


neously, are special cases of our algorithm. Each one of them can be obtained by appropriately 
choosing the scheduling parameters {T^} and In the sequential algorithm, all links 

update their own strategies sequentially, whose scheduling parameters {T^} and are 
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Algorithm 1 Asynchronous Distributed Energy-Efficiency (ADEE) Algorithm 
1 : Set t = 0 and initialize Qfc(O) G Wk^k; 

2: FOR t = 0 : T^ax 
FOR k = l: K 
IF t G 

Update Qfc(t+1) using Algorithm 2 in the following subseetion based on Q_fc(r^(f)). 
ELSE 

Qk{t + 1) = Qk{t), 

END 

END 

END 


ehosen as = {mK + G J\f+}={k,K + k,2K + fc, • ■ ■ = t,\/k,r, where 

N'+ denotes the diserete set Af+ = {0,1, 2, • ■ ■ Similarly, in the simultaneous algorithm, 
all links update their own strategies simultaneously with the seheduling parameters ehosen as 


= t,\/k,r. In eontrast, for the DBEEE [15| and TEEE algorithms [16|, eaeh 
transmitter should eolleet the knowledge of the ehannels to all the reeeivers over the network 
at the initial stage of the algorithm. Moreover, in eaeh iteration all the reeeivers should feed 
baek the eomplex-valued matriees to all the transmitters in the network, whieh is not sealable 
for large seale networks. 


A. Per-link Response Problem 

The best response problem for link k ean be written as 

Cfc(Qfc) 


max . , p . 

QfeSWfc ^"^iQk)+Pc 


(23) 


Sinee the numerator and denominator in the objeetive funetion of Problem (23) are eoneave 


and affine in respeetively, the objeetive funetion of (23) is a pseudo-eoneave funetion p3|. 
Moreover, the eonstraint in ([23]) is eonvex. Henee, problem ([23]) ean be solved by the following 


lemma, the proof of whieh ean be found in Proposition 6 of [411. 


^As Q_fc is treated as constant, Ck{Qk) is used to represent CkiQk, Q-fc) for simplicity. 
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Lemma 4: Define function G{k) as 

G{k) = max log 2 |I + - K(tr(Qfc)+Pc), 


QfcSWfe 


(24) 


For fixed k, the solution of (24) is denoted as QKk). Then, solving (23) is equivalent to finding 
the root of the equation G{k*) = Gk{Q.l{K,*)) — K*(tr(Q^(K*))+Pc)=0. □ 

Lemma 4 gives us insights to solve (|23]). We can solve Problem ([24]) with fixed k firstly, while 


the optimal k can be searched via the Dinkelbach method [42|. 


Now we attempt to solve (24) for fixed k. To this end, we first write the eigenvalue decom¬ 
position (EVD) of for each k G as 


= UtDjUj, 


(25) 


where Ufc is a semi-unitary matrix of the eigenvectors with n = rank(Hf,R-iHfc,fc), Bj, G 


^rkXTk 


is a diagonal matrix with = diagjdfc^, • • • , d/c.r*} being the eigenvalues. 


Then, given k G and Q_fc G >V_fc, the solution to problem ( [24| ) with fixed k is p^ : 

Qk = UfeAfcUf, 

where = diagjgfc^, ■ ■ ■ , gk,rk} represents the power allocations on all subchannels with 

1 1 


(26) 


Qk,m — 




(27) 


(k P Afc)ln2 dk,r, 

where Afc > 0 is the Lagrange multiplier associated with the power constraint, which should be 
chosen to satisfy the complementarity slackness condition: A^ (tr(Qfc) — Pp) = 0. 

After solving problem ( [2^ , we utilize the Dinkelbach method to update n as follows 


f^n+1 — 


tr(Ql''^*(Kn))+Rc’ 

where n is the iteration index. 

To summarize the above analysis, we give the following algorithm to solve the per-link problem 


-,n = 1, 2, • • • , 


(28) 


in (23). 


Algorithm 2 The Dinkelbach method to solve ( 

23 

) 

1: Initialization: kq satisfying G{Ki 

2 : For given solve (24) to get 

3: If (7(«:„) > e, update Kn+i in ( 

))' 

he 

28 

> 0, tolerance e, iteration number n = 0; 
optimal Q^"'^*(k„); 

), n = n + 1, go back to step 2. Otherwise, terminate. 
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Remark 2 - Distributed nature of the algorithm: For transmitter k to update Q^, it requires 
receiver k to feed back the channel matrix ^ and the IPN covariance matrix according 


to (24). It is well-known that the channel matrix can be estimated at the receiver. For 


the IPN covariance matrix R^, it can be easily computed at the receiver as follows: First, 
transmitter k sends a sequence of training sequence to receiver k. The statistical information of 
the training sequence is assumed to be known at the receiver. Then the IPN covariance matrix 
can be obtained by subtracting the covariance matrix of the received training sequence from the 
covariance matrix of the total received signal Hence, the ADEE algorithm can be performed 
in a totally distributed and asynchronous way without the need of information exchange among 
different links. In contrast, for the DBFEE p3| and TEEE algorithms [16|, in each iteration all 
the receivers should feed back the complex-valued matrices to all the transmitters in the network, 
which is not scalable for large scale networks. 


B. Convergence Analysis 

In this section, the sufficient condition that guarantees the global convergence of the ADEE 
algorithm is given. Interestingly, we find that this convergence condition is the same as the 
uniqueness condition obtained in Theorem 2, as proved in the following. 


Theorem 3: In the case of full column-rank channel matrices, suppose that condition (18) in 
Theorem 2 is satisfied and the maximum transmit power is large enough, then as Tmax —>■ oo, 
the sequence generated by the ADEE algorithm converges to the unique NE of game Q, for any 
given set of feasible initial conditions and updating schedules. 

Proof: Please see Appendix □ 

The proof for the more general case can be derived similarly. It is omitted for brevity. 
Remark 3 - Robustness of the algorithm: The condition for the convergence of the ADEE 
algorithm is independent of the update schedule for each link. Hence, all special cases of the 
ADEE algorithm like the well-known sequential and simultaneous updating are guaranteed to 
converge under the same condition in ( [T8] ). Eurthermore, there is no restrict constraints on the 
updating schedule on each link, so that some links are allowed to update their transmit covariance 
matrices more often than others, without affecting the convergence of the algorithm. 


^Note that this has been done in most of the existing distributed algorithms 


0 - 0 - 0 ’ 0 - 0 - @’ 0-0 
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V. Performance Analysis of the Proposed Algorithm 

In this section, we give the performance analysis for the proposed algorithm. Firstly, the impact 
of the circuit power consumption on the system performance is studied. Then, we investigate 
the tradeoff between SE and EE. 

A. Impact of Circuit Power Consumption on the System Performance 

In this part, we investigate the impact of circuit power consumption on the overall SE and EE 
performance. Due to the coupling interference among different links, it is difficult to provide the 
analytical results. To simplify the analysis, we only consider the extreme case that all the links 
are separated far away with each other so that the interference among all the links reduces to 
almost zero. In this case, the overall system can be regarded as K independent links. Then, we 
only need to study the impact of circuit power consumption on the each link fc’s SE and EE 
performance, which is given in the following theorem. 

Theorem 4: When the links are separated far away and no power constraints are imposed at the 
transmitters, the maximum achievable EE of each link k for a given circuit power consumption 
Pc decreases with Pc, but the corresponding SE of each link k increases with Pc- 

Proof: Please see Appendix □ 

To validate the analysis in Theorem 4, we plot a figure in Eig. [^to show the impact of circuit 
power consumption on system’s SE and EE for the same scenario in Eig. with Dcmss = 5. The 
sum-EE and the sum-SE are defined as EEgum = EE*,, SEsum = Yl!k=i ^k, respectively. 

Eor illustration purpose, the y-axis is shown by lOlog^^gEEsum- It can be seen from the figure 
that the sum-EE decreases with the circuit power consumption while the corresponding sum- 
SE increases with it for all considered numbers of antennas, which validates the correction of 
the theorem. Erom the simulation section, the above trend also holds for the general case with 
sizeable interference among the links, though the above theorem is derived for the extreme case. 

B. Tradeoff between Spectral and Energy Efficiency 

In this part, we study the tradeoff between SE and EE for the proposed ADEE algorithm and 
the asynchronous distributed spectral efficient (ADSE) algorithm pT| . Generally, it is difficult to 
study this problem analytically due to the coupling interference. Eor this reason, we investigate the 
tradeoff for two special cases: low transmit power constraint and high transmit power constraint. 
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Fig. 2. The impact of circuit power consumption on system’s SE and EE when links are separated far away. 


1) Low transmit power constraint: Pt —^ 0.’ In this case, the interlink interferenee ean be 
negleeted, i.e., Rfc ~ I, Vfc, and the network reduees to K independent point-to-point links. 
Henee, we only need to analyze one link’s performanee. Without loss of generality, we only 
foeus on the performanee of link k. Denote the transmit power of link k of the ADEE algorithm 


as Pfc G [0, Pt]. Sinee Pp —)■ 0, from (E.3) there is only one stream for link k that transmits with 
positive power. Henee, the SE for link k ean be written as Ck (Pk) = log 2 (1 + dk,iPk), where 
dk,i is the maximum eigenvalue of ii^fjik,k- The optimal EE for link k ean be written as 

Ck {Pk) 


EEfc (Pk) — 


Pk + Pc 


(29) 


Then, we have 

where EE'^ (P^) is the derivative of EE*. (P*.) w.r.t. Pk- Henee, EE^ (P^) is inereasing for Pk G 
[0,Pt]. Then, the optimal transmit power is P^ = Pp, i.e., link k transmits with its maximum 
power. 

Eor the ADSE algorithm, eaeh link maximizes its SE selfishly by always using its maximum 
power, i.e., P^ = Pp. Henee, both the ADEE algorithm and ADSE algorithm aehieve the same 
SE and EE. 
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2) High transmit power constraint: Pt oo; In this case, the interlink interferenee eannot 
be negleeted any more. Due to the interferenee, the analysis beeomes more diffieult. To simplify 


the analysis and get insights, we eonsider a symmetrie system similar to one in |j^| for the SISO 
frequeney seleetive interferenee ehannel. We assume that H = Vfc and ^ 
where a is a eonstant. In this symmetrie system, all links transmit with the same eovarianee 
matriees, i.e., Q = Qfc , 'ik. Denote the transmit power as tr(Q) = P. 

The overall network EE is 

Plog2 ll + H^R-^HQl 


EE 


P + Pc 


and the overall network SE is 


SE =Plog211 + H^R^'HQI , 


(31) 


(32) 


where R = I + (P - l)a2HQH^. 

Eor the ADSE algorithm, eaeh link selfishly maximizes its own SE by using its maximum 
transmit power, i.e., tr(Q) = P = Pp. When Pp —)■ oo, the interlink interferenee beeomes 


very large. Aeeording to Theorem 3 of [241, one of the optimal solutions employs beamforming 
(1-stream signaling) for all links. Henee, the optimal eovarianee matrix ean be written as Q = 
Ppww^, where w is the beam direetion with unit norm. Then, the network SE of the ADSE 
algorithm in the high transmit power is upper bounded by 

5'-Pt \ f. 1 


SEadse = lim A:iog2 1+ . , p 

Pt^oo y 1 -|- — Ij^fPp 


= ^log2 1 + 


a^{K-l) 


(33) 


where g = w^H^Hw is a constant. The corresponding EE is EEadse = lim SEadse/(-Pt + -Pc) 

Pt^oo 

0, whieh is not desirable from the EE point of view. 

On the other hand, for the ADEE algorithm, the system will not use the maximum power to 
transmit due to the power value in the denominator of pTj ). Henee, the interlink interferenee may 
not be so large. Thus, the above derivations for the high interferenee seenario are not applieable 
and it is diffieult to analyze the performanee of the ADEE algorithm in this ease. However, 
from the above diseussion, we ean eonelude that its overall EE will inerease when the transmit 
power eonstraint is low, and keeps eonstant when the transmit power eonstraint is high, whieh 
is signifieantly larger than that of the ADSE algorithm. The overall SE has the same trend as 
the overall EE. 
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VI. Simulation Results 

We consider an ad hoc network contained in a 250m x 250m square area, in which all links are 
randomly distributed. The distance from one transmitter to its unintended receiver is at least 35m. 


The channel is modeled by path-loss [431 and independent Rayleigh fading with the complex 


normal distribution, CA/'(0,1). The channel path-loss is modeled as 38.46 -f 351ogio(d) [431. 
Each channel realization is obtained by generating a random set of links’ positions as well as 
fading channel realizations. It is assumed that the transmitters and the receivers have the same 
number of antennas. Unless otherwise specified, the other main system parameters are given in 
Table For comparisons, the metrics used are the sum-EE and the sum-SE, which are defined 
as EEgmjj I —1 EE/j, SEgm^ ^^k=i respectively. 


TABLE I 

Main simulation parameters 


Parameters 

Value 

Number of links K 

4 

Number of transmit antennas M 

4 

Number of receiver antennas N 

4 

Direct-channel distance Udirect 

80 m 

Noise power 

-106 dBm 

Circuit power consumption Pc 

23 dBm 44 

Tolerance e 

10 (- 5 ) 

Maximum transmit power Pt 

30 dBm 

Maximum number of iterations Tmax 

20 

Channel path loss model 

38.46-f351ogio(d) 43 


A. Convergence behavior of the ADEE algorithm 

Figj^ illustrates the convergence behavior of the ADEE algorithm for different updating 
schedules for one randomly generated channel realization. For comparison, the performance 
of another updating scheme, named ‘unbalanced ADEE’, is also shown, where the parameters 
are chosen as = {k, 2k, 3k, • • ■ T^(f) = t, \/k, r. Hence, in the unbalanced ADEE, Link 1 

is set to update its strategy faster than the other three links, and Link 4 is set to be the slowest 
to update its strategy. Also, in each time, the number of links updating their strategies is no 
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Fig. 3. Converegence behaviors for different updating schedules. 



Fig. 4. Average sum-EE (left) and sum-SE (right) versus Pt for various algorithms. 
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less than that of the sequential ADEE (only one link), no more than that of the simultaneous 
ADEE (all four links). To make the figure not too erowded, we report only the eurves of two 
links (Eink 1 and Eink 4). It can be seen from this figure that the simultaneous ADEE converges 
faster than the other two schemes, and can converge within a few iterations. However, it takes 
about 12 iterations for the sequential ADEE to converge. The reason is that each user in the 
sequential ADEE is forced to wait for all the users scheduled in advance, before updating its 
own strategy. Moreover, the unbalanced ADEE converges a little faster than the sequential ADEE 
due to more links involved in updating their strategies in each iteration. Erom this figure, we 
find that different scheduling methods yield almost the same performance. Eor this reason, we 
only report the performance of the simultaneous ADEE due to its rapid convergence speed in 
the following simulations. 


B. Performance Comparison with Existing Algorithms 

We next compare the performance of the ADEE algorithm with some existing algorithms, 
including the ADSE algorithm pT| where each link always uses its maximum transmit power, 
and DBEEE algorithm p3| . Since in the proposed ADEE algorithm each link only attempts 
to maximize its own EE selfishly, its achieved overall EE is generally suboptimal. Hence, it is 
interesting to study the performance gap between the proposed ADEE algorithm and the (near- 
)optimal sum-EE maximization algorithm. Eor this reason, we simulate the TEEE algorithm 


in [161 that aims at the sum-EE maximization. Moreover, the WMMSE algorithm in [14| that 
focuses on the sum-SE maximization is also simulated. 

1) Impact of the Maximum Transmit Power on the System Performance: Eig. shows the 
average sum-EE (left) and the average sum-SE (right) versus the maximum transmit power for 
various algorithms. It can be seen that in the low transmit power regime, i.e., Pt < 10 dBm, 
both the sum-EE and the sum-SE of all algorithms monotonically increase with the maximum 
transmit power. Interestingly, we find that in this regime the proposed ADEE algorithm almost 
achieves the optimal sum-EE and sum-SE. That is, the performance gap with the optimal sum-EE 
and sum-SE performance approaches zero. Eortunately, many short-distance wireless networks 
operate in this regime, such as ad hoc networks [[43|, [|46|, femtocells, and wireless sensor 


networks [47|. However, in the high transmit power regime, the sum-EE of the ADRM algorithm 


decreases dramatically, while its sum-SE performance becomes stable. This is because in the 
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high transmit power regime, the system beeomes interferenee limited and inereasing the transmit 
power may slightly help inerease the sum-SE, whieh leads to a signifieant performanee loss in 
terms of the sum-EE due to the inereased transmit power. Note that the sum-SE aehieved by the 
ADEE algorithm is eomparable with that aehieved by the ADRM algorithm, whieh eorroborates 
with the analysis in Seetion |Vj Eig. also shows that the ADEE algorithm outperforms the 
DBE algorithm in terms of both the sum-EE and the sum-SE. This is beeause the DBEEE 
algorithm is primarily designed for the symmetry system and may not be suitable for the 
asymmetrie ad hoe network eonsidered here. Moreover, this algorithm requires a substantial 
feedbaek overhead on the network and further the power eonsumption due to these information 
exehanges is not eonsidered. As expeeted, both the sum-EE and the sum-SE performanee of 
the TEEE algorithm is superior to that of the ADEE algorithm in the high transmit power 
regime due to the selfish nature of the ADEE algorithm. However, this benefit eomes at the eost 
of the heavy information exehange overhead, high eomputational eomplexity, synehronization 
requirements of the networks, and the need for eoordination among the users. In eontrast, the 
proposed ADEE algorithm does not need information exehange among different links, and ean 
be implemented in a totally distributed and asynehronous manner, whieh is appealing for some 
praetieal applieations. 

2) Impact of the Number of Antennas on the System Performance: Eig. illustrates the 
average sum-EE (left) and the average sum-SE (right) versus the number of antennas for various 
algorithms. As expeeted, inereasing the number of antennas at both the transmitters and the 
reeeivers leads to a signifieant inerease in both the sum-EE and the sum-SE due to the faet that 
a larger number of degrees of freedom (DoE) in the spatial domain ean be exploited to strengthen 
the signal power reeeived by the intended reeeiver while avoiding the interferenee imposed on 
the unintended ones. It ean be seen from this figure that the ADEE algorithm outperforms the 
DBE algorithm in terms of the sum-EE and the sum-SE. The performanee gain monotonieally 
inereases with the number of antennas, meaning the ADEE algorithm uses the spatial resourees 
more effeetively. Again, the ADEE algorithm is observed to aehieve a signifieantly higher sum-EE 
than that aehieved by the ADRM algorithm and the performanee gain beeomes more signifieant 
with more antennas. The sum-SE of the ADRM algorithm is slightly higher than that of the 
ADEE algorithm. As expeeted, the TEEE algorithm has superior performanee over the ADEE 
algorithm in terms of the sum-EE and sum-SE. 
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Number of Antennas Number of Antennas 


Fig. 5. Average sum-EE (left) and sum-SE (right) versus the number of antennas for various algorithms. 



Fig. 6. Average sum-EE (left) and sum-SE (right) versus the direct-channel distance Ddirect for various algorithms. 
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3) Impact of the Direct-Channel Distance on the System Performance: Fig. [^illustrates the 
impaet of inereasing the direet-ehannel distanee /indirect on the sum-EE and sum-SE. It ean be 
seen from this figure that both the sum-EE and sum-SE monotonieally deereases upon inereasing 
^direct- The reason is that increasing /indirect will degrade the direct-channel gains, and thus 
reducing the attainable sum-EE and sum-SE. This figure shows much of the same properties as 
Eig. Eor example, the ADEE algorithm outperforms the DBEEE and ADRM algorithms in 
terms of the sum-EE, and the sum-SE loss is negligible compared with the ADRM algorithm. 
With much more information exchange overhead and higher computational complexity, the TEEE 
and WMMSE algorithms have the best sum-EE and sum-SE performance, respectively. 

4) Impact of the Circuit Power on the System Performance: The impact of the circuit power 
on the attainable sum-EE and sum-SE is characterized in Eig. [^ Eor illustration purpose, the 
y-axis is shown by lOlogj^oEEsum- As expected, the sum-EE achieved by all the algorithms 
decreases with the increase of the circuit power due to the increased power consumption and 
finally converges to almost the same value. On the other hand, the corresponding sum-SE of the 
ADEE algorithm and the DBEEE algorithm monotonieally increases with the increase of circuit 
power and finally converges to the sum-rate achieved by the ADRM algorithm. This indicates 
that more circuit power will encourage each link to use more power to achieve a higher SE, 
which shows the trend in Theorem 4 also holds in the more general case. Note that the similar 


trend has been observed for the MIMO interference channel in [ 151. It is interesting to find that 
the sum-EE gain provided by the TEEE algorithm over the ADEE algorithm is small and fixed 
during the overall circuit power consumption regime. 

5) Impact of number of links on the System Performance: Einally in Eig. we investigate the 
impact of increasing the number of links on the sum-EE and sum-SE performance of different 
algorithms. It can be seen from this figure that the sum-EE of the ADEE algorithm decreases with 
the increase of the number of the links. The reason is that when the number of links increases, the 
interference power received at the receivers increases. To compensate for this negative effects, 
each transmitter will increase its transmit power. However, this slightly improves the sum-SE 
but significantly reduces the sum-EE. To improve the performance of our algorithm for this 
interference-limited scenario, we may combine our algorithm with some scheduling methods, 
which is beyond the scope of this paper. On the other hand, by increasing the number of links, 
the sum-EE performance of the DBEEE algorithm gradually increases. This can be explained as 
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Fig. 7. Average sum-EE (left) and sum-SE (right) versus the circuit power for various algorithms. 
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Fig. 8. Average sum-EE (left) and sum-SE (right) versus the number of links for various algorithms. 
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follows: As the number of links increases, the interference power perceived at different receivers 
become comparable with each other due to the law of large numbers. As a result, the system 
will become symmetry, for which the DBFEE algorithm will be suitable. It can also be observed 
from this figure that there is a point beyond which the DBEEE algorithm outperforms the ADEE 
algorithm. However, this superiority comes at the cost of more energy consumption incurred by 
the information exchange, which is not accounted for in the simulations. Einally, it is observed 
that both the sum-EE and sum-SE increases with the number of links due to the multiuser 
diversity. 


VIE Conclusion 

In this paper, we have considered a game theoretical formulation of the maximization of 
the EE on each link, subject to the power constraints, in the MIMO interference channel. 
We have provided a complete characterization of the game, by showing the existence of the 
NE and deriving the sufficient conditions for the uniqueness of the NE for the case of large 
enough maximum transmit power constraint. We have provided a totally asynchronous iterative 
distributed algorithm, named the ADEE algorithm, to reach the NE of this game. The ADEE 
algorithm has three advantages: no rigid scheduling in the updates of the players is required, the 
synchronization requirement of both sequential and simultaneous algorithms can be removed and 
information exchange among the players is not necessary, which make the proposed algorithm 
appealing to implement in practice. Furthermore, the sufficient conditions guaranteeing the 
convergence of the algorithm have been provided. Interestingly, our sufficient condition does 
not depend on the scheduling of the links. Extensive simulation results have shown that the 
ADEE algorithm performs better than the DBFEE and ADRM algorithms in terms of the sum-EE 
performance, and is comparable to the ADRM algorithm in terms of the sum-SE performance. 


Appendix A 
Proof of Eemma 1 


Before proving Lemma 1, the following definition is given. 


Definition 2 [48 A real-valued differentiable function /(X) : X 
quasiconcave on if its sublevel sets Sa = {X G §:!J^^^|/(X) 

(strictly) convex. 


G is (strictly) 

> a}, for a G R, are 
□ 
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For simplicity, we omit the dependeney of EEfc(Qfc, Q_fc) and Ck{Qk, Q-k) on Q_fc. Define 

two funetions /i(Qfc) = Ck{Qk) and /2(Qfc) = tr(Qfc)+Pc, and define /(Qfc) = EEfc(Qfc) = 

/i(Qfc) 

f2{Qk)- 

Now, we give the proof of the lemma. Sinee / 2 (Qfc) > 0, Sa is equivalent to 


= {Qk e - «/2(Qfc) > 0}. (A.l) 

Sinee f 2 iQ,k) is affine in Q^, Sa is (strietly) eonvex for any given a if fi{Q,k) is (strietly) 
eoneave in Q^. In the following, we will prove /i(Qfc) is eoneave w.r.t. for any ehannel 
matrix ^ and is strietly eoneave w.r.t. if ^ is full eolumn rank. 


To this end, we adopt the teehnique in [24|. We eonsider the eonvex eombination of two 
different feasible G Wfc and G Wfc, whieh is 


Qfc = fZfc + (1 — f)Xfc = Xfc + fYfc, (A. 2) 

where 0 < f < 1, and — X^, whieh is a non-zero Hermitian matrix. Obviously, 

Qfc G Wfc. Then /i(Qfc) is (strietly) eoneave w.r.t. Qfc if and only if (d^/df^)/i(Qfc) < 0 
((dVdf2)/i(Q,) < 0) for any feasible Xfc and Zfc. Denote Afc = H|^^R^^Hfc fc, the seeond 
derivative of /i(Qfc) w.r.t. t can be ealeulated as 

= -j^tr (AfcYfc(I + AfcQfc)-'AfcYfc(I + AfcQfc)"') . (A.3) 

Sinee (I + AfcQfc)”^ G it is possible to write (I + AfcQfc)”^ = DfcD^ with G 

Then, it is seen that 

= -j^tr ((Df AfcYfcDfc) (Df AfcYfcDfc)) < 0. (A.4) 

Furthermore, when Hfc fc is a full eolumn-rank matrix, AfcYfc is always non-zero sinee Yfc is a 
non-zero matrix. Then (d^/df^)/i(Qfc) < 0 always holds. Henee, the proof is eompleted. 


Appendix B 
Proof of Theorem 1 

We prove the existenee of NE for game Q using the following well-known game theory result. 
Lemma 5 [49J : The strategie noneooperative game Q = {^fcjfce^^^> {^fclfce^xl admits 
at least one NE if, for all k G 'F/c:!) the set Wfc is a nonempty eompaet eonvex subset of a 
Euelidean spaee, and 2) the payoff funetion Uk is quasieoneave on Wfc. □ 

It ean be easily verified that game Q satisfies the two eonditions required by Eemma 5 
aeeording to Eemma 1 and the eonvexity of the admissible power set. 









Appendix C 
Proof of Theorem 2 
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Before proving Theorem 2, we provide three lemmas [38| that will be used in our derivations. 
Lemma 6: Given matriees A G G the relation ||AB ||2 < ||A|| 2 ||B ||2 holds. □ 

Lemma 7; If A,B G then we have p(AB) < p(A)p(B). □ 

Lemma 8: Given matrix A G with eigenvalues {Xi,i = 1, • • • ,n} and eigenveetors 

{xj, i = 1, • • • , n}, and matrix B G with eigenvalues {/x*, i = 1, ■ ■ • , n} and eigenveetors 

{yj, i = 1, • • • , n}, the eigenvalues of A ® B are given by {Aj/ij, x = 1, • • • , n; j = 1, • ■ ■ , m} 

and the eorresponding eigenveetors are given by Xj (g) y^-, x = 1, • • ■ , n; j = 1, • • • , m. □ 

Based on the above results, we now proeeed to prove Theorem 2. Our main idea lies in 


showing that under eondition (18), eondition (12) in Lemma 3 will hold. To this end, we will 


employ the mean-value theorem for eomplex matrix-valued funetions stated in [311. That is, for 
any two different points and Q^]., there exists some t G (0,1) sueh that 




( 2 ). 


<|Dq_.F»(A) 


qA- 


■q3 


(C.l) 


where A = -I- (1 — f)Q'Lfc and Dq ,^Ffc(A) represents the Jaeobian matrix of funetion 


v(2) 


Ffc(Q-fc) w.r.t. Q-fc, evaluated at Q_fc = A [[39|. 

Let DR,j^Ffc(A) be the Jaeobian matrix of funetion Ffc(Q_fc) w.r.t. Rfc, evaluated at = 
Rfc(A). Let Dq ,^Rfe(A) be the Jaeobian matrix of funetion HkiQ-k) w.r.t. Q_fc, evaluated at 
Q_fc = A. Then, we have: 


Ft(QL‘() - F,(Q«) 


( 2)2 


< Dr, Fi.('A) ■ Dr 


Rfe(A) 


QX - Q 


( 2 ) 

-k 


<||Dr,F,(A)||2||Dq_,R,(A) 


qA 


qS 


(C.2) 

(C3) 


where ( |C.2| ) follows from the ehain rule and ( |C.3| ) follows from Lemma 6. 

Our goal is to obtain the upper bound of ||DRj^Ffc(A)|| 2 ||DQ ,^Rfc(A)|| 2 . If this upper bound 
is less than ^yi/{K — 1), then Condition (12) is satisfied. In the following, we derive the upper 
bounds of ||DRj^Ffc(A )||2 and ||DQ_^Rfc(A) 11 respeetively. 


A. The upper bound of ||Dq j^Rfc(A)||2 

To derive the upper bound of ||Dq ,^Rfe(A) we should obtain the expression of Dq ,^Rfc(A). 
To this end, the expression of the Jaeobian matrix Dq ^.R^, whieh is a funetion of Q_fc, should 
be obtained firstly. Then Dq ,.Rfc(A) ean be obtained by inserting Q_k = A into DQ_j,Rfc. 
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The Jacobian matrix DQ_^Rfc can be computed using the three-step procedure in [39|. Specif¬ 


ically, we compute firstly the differential of and then the Jacobian matrix. Function R^ is 
differentiable at Q_fc, with differential given by dR^ = By vectorizing dR^ 

and using the equality vec(ABC) = (C^ ® A)vec(B) Bm, we obtain 


dvecRfc=[Hi , 










K,k 


iiK,k] dvec(Q_fc) 


(C.4) 


Using the identification rule in |39|, we can obtain Dq ,.Rfc given in ( [T7| ) in Theorem 2. Note 
that Dq ,^Rfc does not depend on Q_fc. Hence, we have Dq ,^Rfc(A) = Dq ,^Rfc. 

B. The upper bound of ||DRj^Ffc(A )||2 

Similarly, to obtain the upper bound of ||DR,j^Ffc(A)|| 2 , we should obtain the expression of 
DRj.Ffc, which is a function of Q_fe. Then the expression of DRj,Ffc(A) can be obtained by 
inserting Q_fc = A into it. 

According to Lemma 1, the EE function EEfc(Qfc) is strictly quasiconcave on since 

is assumed to be full column rank. Also, the maximum transmit power is assumed to be 
very large. Then, for given Q_fc (and thus R^), the gradient of EEfc(Qfc) w.r.t. Q^, evaluated 


at Qfc=Ffc, must be zero [|^. Thus, using d |BXC|/aX = |BXC| [C(BXC)^^B] [jsy, we 
have 

G(Rfc,Ffc)^Hj,(R, + Hfc,fcFfcH^,)-iHfc,fe-«J = 0. (C.5) 

where ak = Ck{Q,k,Q,-k)/{Pk + Pc)- Since function EEfc(Qfc) is strictly quasiconcave in Qfc 


for given R^, there exists a unique globally optimal solution F^ that satisfies (C.5) [501. Hence, 


equation (C.5) defines an implicit function [52|. Taking the derivative of (C.5) w.r.t. R^ and 
using the chain rule, we have 


Dr. G -f Dp, G ■ Dr. Ffc — 0. 


(C.6) 


Now we first obtain the Jacobian matrices Dr^^G and Dp^^G, then Dr^^F/j can be solved from 


(C.6). We again use the three-step procedure in [39| to compute Dr^^G and Dp^^G. Function 


G is differentiable w.r.t. both R^ and F^, with the differential given by 


dG — —J;i;dRfc 




(C.7) 


''For simplicity, the dependency of function G(Rfc,Ffe) on both and Fa, is omitted. 
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where = (R^ + and we used dX ^ = —X MXX ^ |39|. Then by vee- 

torizing both sides of (|C.7|) and using the identifieation rule in [|39|, we ean obtain the Ja- 


eobian matriees Dr, G = — Ufc and Df.G = — , with (8) and 


IH 


k 


^k,k^k 


Vfc = 


I fc. Note that is a Hermitian matrix. Sinee the ehannel matrix 


Hfc fc is assumed to be full eolumn rank and Jfc is nonsingular, ^ is nonsingular. Then 


Vfc is also nonsingular. This is beeause [ [5T | 

(A0B)^^ = , for all nonsingular A, B. 


(C.8) 


Henee, Dr F^. ean be solved from (C.6): 


Dr, F, = -vr^u,. 


(C.9) 


Thus, DRj^Ffc(A) ean be obtained by inserting Q_fc = A into (C.9): 


Dr.F»(A) = -V,-‘(A)Ut(A 


(C.IO) 


Then, we have 

||Dr,F,(A)||^ = p(Uf (A)V,-2(A)U,(A)) < p(V,-2(A))p(U.(A)Uf (A)), (C.ll) 

where the inequality follows from Lemma 7. We now aim to obtain the upper bound of p(V^^(A)) 
and p(Ufc(A)U|^(A)), respeetively. 

1) The upper bound of p(V^^(A)).- The upper bound of p(V^^(A)) ean be obtained as 
follows: 


P(V,-^(A)) 

= P ((H^tJ^AjH^ 0 H"jJ,(A)Ht,i)-") 


< 


, -4 


|A„,„ + H 


-4 


2 = 1 


(C.12) 

(C.13) 

(C.14) 

(C.15) 

(C.16) 

(C.17) 


where Jfc(A) is obtained by inserting Q_fc = A into J^, (C.14) follows from (C.8) and the 


equality (A (g) B) (C ® D) = AC ® BD [511, (C.15) follows from Lemma 8, (C.16) results 
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from the fact that H|^Jfc(A)Hfc G the last inequality results from Amin(A^BA) > 

Amin(A^CA) for all B ^ C ^ 0 and the following relations 


= I + E ah" + Ht,tF,(A)H"^ 


i^k 




I + tr (Ffc(A)) tr (A) 


H 

i,k 


i^k 


K 


-1 


[ffc 


(C.18) 

(C.19) 

(C.20) 


E I + Pr5^H,,fcH 

V i=i 

2) The upper bound o/p(Ufc(A)U^(A)).- The upper bound of p(Ufc(A)U|^(A)) can be 
obtained as follows: 

p(Ui(A)Uf (A)) = p ((HJ^JPA)) 0 H"jJt(A)))(J:(A))H;,, 0 jf (A))Ht.»))(C.21) 
= p(H"iJ"(A)Ji(A)HJ,i0H",J,(A)J«(A)H,,J (C.22) 

= (p(H"tJj(A)Hpt))" (C.23) 

<(p(H"tHtp.))" (C.24) 


where in (C.22) we use (A 0 B) (C ® D) = AC 0 BD [511 , (C.23) follows from Lemma 8, 


and (C.24) follows from J^fA) P I. 


Combining ( jC.llj ) with ( jC.lVj ) and ( |C.24| ), we obtain 




||DR,Ffc(A)|L<- 
Finally, by combining ( |C.3| ) with ( [T7| ) and ( |C.25| ), we have 


Fi(QL‘l) - F^(Q«) 


(2)^ 


< ttfc 


QA - Q 


(C.25) 


(C.26) 


where au is given in ( [T^ . Hence, if condition (18) is satisfied, condition ( [T^ in Lemma 3 holds, 
which completes the proof. 


Appendix D 
Proof of Theorem 3 

Before proving the theorem, we introduce some basic definitions that will be used in our 
derivations. Given the multiuser mapping F(Q) defined in (j^, we introduce the following block- 
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maximum norm on C 


KMxKM 


, defined as [53| 


l|F(Q)||i7’block ~ ||Ffc(Q-fc)||i?- 


Let IN 


and let 


be the veetor maximum norm, defined as [38| 

Halloo,vec = |a^fc|,X e I 




(D.l) 


(D.2) 


I oo mat denote the matrix norm indueed by 

K 

IIAI 


«,vec’ given by 1^: 

loo,mat - max^|[A];,^J,A e 
r=l 

Finally, we introduee the nonnegative matrix B G defined as 

al, if fc 7^ r, 


(D.3) 




0, otherwise. 


(D.4) 


Define Omax = max al and /3 = {K — l)amax < 1, then ||B||^ ean be easily eomputed as 
fcei'K 


|B| 


oo.mat 


= / 3 . 


(D.5) 


Based on the above results, we then give the eontraetion property of the multiuser mapping 
funetion in the following lemma, whieh will be used in the proof of the theorem. 

Lemma 9: If the uniqueness eondition in ( [T^ in Theorem 2 is satisfied, then the multiuser 
mapping funetion F(Q) defined in (|^ is Lipsehitz eontinuous on W: 


F(Q(^)) -F(Q 


(2)^ 


IF,block 




( 2)1 


1F,block 


(D.6) 


VQ*^L^q( 2) g where ||-Unblock is defined in (|D.1|). Furthermore, the mapping is a bloek- 
eontraetion with modulus \/^ < 1- 

Proof: Given ■ ■ ■ ,Q^^) G W and = (Qf\ ■ ,Q^^) G W, define, for 


eaeh k G = 


F^(Qtl) - F»(Q™) 


{2)^ 


and Cfc = 


Qf’ - Qf 


. Then we have 


eF.= 


Ft(QL'))-Fi(Q'4) 


(2)^ 






Qf’-Qf >EP]c/i(D-7) 






VQ(L,Q(2) e yy and Wk G ^ k - 

Introdueing the veetors ep = ■ ■ ■ , and e = [ei, • • • , the set of inequalities 


in (D.7) ean be rewritten as: 


0 < Oi. < G W. 


(D.8) 
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Then we have |53| 


lepll < llBell < ||B|| Jlell = /3||e| 

I -^lloo,vec — II lloo,vec — II lloo,matll lloo,vec <^11 I 


(D.9) 


where ||||^_^^, and 


I oo.mat 


are defined in (D.2) and (D.3), respeetively, and the last equality 


follows from (D.5). 


Finally, using (D.9) and (D.l), one obtains 

2 


||F(Q<»)-F(Q<"l)||;^ 

G W. Henee, the lemma follows. 


block ll®^lloo,vec 




( 2)1 


IF, block 


(D.IO) 


□ 


Interestingly, we find that if eondition (18) in Theorem 2 is satisfied, the multiuser mapping 
function F(Q) is not only a block-contraction for the Frobenius norm as shown in Lemma 2 
but also block-maximum norm as shown in Lemma 9. 


The remaining task is to show that, under condition (18) in Theorem 2, conditions of the 


asynchronous convergence theorem in Prop.2.1 of [401 on page 431 are satisfied. Based on the 


above lemma, the proof can be done in a similar fashion to the proof for the rate maximization 


game in [311 . For this reason, the details are omitted for simplicity. 


Appendix E 
Proof of Theorem 4 

Since the links are located far away, the interference among the links is negligible so that the 
covariance matrix of each link k (i.e., R^) can be regarded as an identity matrix I. Hence, the 
SE function of links k reduces to 


Ck{Pk) = max log 2 I -f . 

tr(Qfe)=Pfc 


(E.l) 


Obviously, the optimal solution of is = UfcPfcU|^, where Ufc is the eigenvector matrix 
of and is nonnegative matrix with P^ = diag{pA:,i, • • • which is the power 

allocation over the streams. By inserting into ( [E.l[ ), the SE function can be equivalently 
obtained as 


rk 


Ck{Pk) = max XI log2(l + 4,mPfc,m) 

m=l 

rk 

Pk,m — Pk-) Pk^m ^ 0 * 

m=l 


(E.2) 
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The solution to the above problem ean be easily solved as = (/^fc — ^ 2, • • ■ , r^}, 

r-k 

where /Xfc is the water-level ehosen to satisfy Pk,m = Pk, and dfc.m’s are the eigenvalues of 

m=l 

whieh are arranged in deereasing order. 

i 

Define gi = Id^] — d^m and gvk+i = oo- Obviously, gi = 0. Then, the SE function can 

m=l 

be rewritten as 

Ck (Pk) = (^(Pk + J2 e [gi, pz+i], / = 1, • • • , r,. (E.3) 

m=l \ \ m=l ) ) 

Note that I represents the number of streams that are allocated with positive power allocation. 

Based on the above analysis, we now provide the properties of the SE function in the following 
lemma. 

Lemma 10: The SE function Ck{Pk) is a continuous, strictly increasing, differential and strictly 
concave function of Pk- Moreover, the first derivative of function Ck{Pk) (denoted as Ck\Pk)) 
is a continuous function of Pk. Also, the SE function Ck{Pk) is twice differentiable in each 
interval of the subregions [p;, , Z = 1, ■ ■ ■ , r^, which is given by 

P'kjPk) = < ^, Pk ^ {gi, gi+i)J = ^, ■ ■ ■ yfk. (E.4) 

Proof: Obviously, in each interval of the subregions, the SE function Ck{Pk) is a continuous, 
strictly increasing, differential and concave function of Pk- We only need to verify the boundary 
points. Eor the boundary points, we have 


lim Cfc(Pfc) = ^log 2 


m=l 

lim Ck{Pk) = 


^k,m 




= lim Ck{Pk), 


Pk^gfi 


= lim Ck{Pk) > 0. 


Pk—^gfi 


(E.5) 

(E.6) 


Equation (E.5) means that the SE function Ck{Pk) is a continuous function and equation (E.6) 
means that the SE function is a strictly increasing and differential function. 

Unfortunately, the SE function is not twice differentiable since 

r riiikTD\ 1 dli_^_i 1 

hm_ Ck [Pk) = -r^^— ^ -r^TTT 
Pk^gi+i 111 2 I ln 2 / + l 


= lim Ck"{Pk). 

Pk^gfi 


(E.7) 


Hence, the concavity of the SE function cannot be proved by using the above method. To 
overcome this issue, we use the duality theory to prove the concavity of the SE function. 


Specifically, for a given Pk, the SE function Ck[Pk) can be obtained by solving problem (E.2), 
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which is a strictly convex optimization problem. It ean also be verified that the eonstraints in 


problem (E.2) satisfy the Slater’s eondition. Henee, the duality gap for this problem is zero. 
Thus, the SE funetion ean also be obtained by solving the following min-max problem 


Tk rk 

Ck{Pk) = min max log 2 (l + 4,mPfc,m) - 

M>0 Pk,m>0,Vm ^^' 

m=l 


rk 

min > 

fi>0 ^ 


log: 


m=l 

dk.m 


m=l 

rk r 

logs 


/xln 2 


E 

m=l 


d, 


k,m 


/i(Pfc)ln2 


rk 

E 

m=l 

rk 

E 

m=l 


1 

1 

h^ 


dk,m 

KPk) 

dk.m 


+ f^Pk 


(E.8) 


(E.9) 


+ fi{Pk)Pk, (E.IO) 


where /r(Pfc) > 0 denotes the optimal value for given Pk, whieh is unique due to the faet that 
problem (E.3|) is a strietly eonvex problem. 


Eor any two points pj^^\plf \ define P^^'^ = OP^^ + (1 — 0)P^‘^’, where 6 G (0,1). Eet 
^,{P^^) and be the optimal /x for Ck{P^'’), Ck{P^k^) and CkiPjf^), respeetively. 

Then, for x = 1, 2, we have 

rk r 


>( 2 ) 


Ck(pi) = E 

m=l 

rk 


logs 


< 


E 

m=l 


log^ 


(- 


dk,m A 

+ rk 

1 

MPi)] 

/i(P^)ln2y_ 

m=l 

In 2 

dk,m 

dk.,m A 

+ rk 

1 


l.i(Pi)\n2j_ 

m=l 

In 2 

dk^m 


+k{p;)pi (E.11) 


+ li(Pl)P't, (E.12) 


where the striet inequality follows due to the faets that fJ^iPk) is not the optimal solution to 


problem (E.9) for given Pl,i = 1, 2, and problem (E.2) has a unique global solution sinee it is 
a strietly eonvex optimization problem [ [^ . By using the above inequality, we have 

(E.13) 


eCk{p^) + (1 - e)Ck{p^,) < Ck{Pi) = Ck{epl + (i - e)Pi). 


Henee, the SE funetion is also a strietly eoneave funetion of Pk. Obviously, the SE funetion 
Ck{Pk) is twiee differentiable in eaeh interval of the subregions [gi, gi^i] , Z = 1, ■ ■ ■ , and the 
seeond-order derivative of Ck{Pk) w.r.t. Pk ean be easily ealeulated in (E.4). □ 

We now proeeed to prove the first part of Theorem 4: the eorresponding SE inereases with 
the eireuit power eonsumption. The EE funetion ean be rewritten as 

Ck{Pk) 


EEi. — 


(E.14) 


Pk + Pc 

Aeeording to Eemma 10, Ck{Pk) is a strietly eoneave funetion of Pk- Moreover, the denominator 


of the EE funetion is affine in Pk. Henee, EE^ is strietly quasieoneave in Pk [481. Then, the 
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optimal solution of Pk (denoted as PI) to maximize the EE funetion is unique [ |54| and should 
satisfy the first order eondition [481: Ck{P^){P^ + Pc) — Ck{Pk) = 0, whieh is equivalent to 

Ck{Pt)-PtC'k{Pt) 


= P 


c- 


(E.15) 


C'k{Pt) 

Our task is to analyze the monotonieity of P^ w.r.t. Pc- If P^ is strietly monotonieally inereasing 
w.r.t. Pc, the eorresponding SE is monotonieally inereasing w.r.t. Pq as well sinee the SE 
is inereasing in Pk based on Eemma 10. However, direetly proving the monotonieity of P^ 
w.r.t. Pc is diffieult. Instead, if the following two eonditions hold: 1) Pc is a eontinuous 
funetion of P^; 2) there is a one-to-one mapping relation between P^ and Pc, this proof ean 


be equivalently transformed into the opposite side Pc is strietly monotonieally inereasing 
w.r.t. P^. Obviously, Pc is a eontinuous funetion of P^ sinee both Cfc(P^) and Ck{Pt) are 
eontinuous funetion of P^ aeeording to Eemma 10. The one-to-one mapping relation between 


P^ and Pc is also obvious: Eor given Pc, there is a unique P^ satisfying (E.15) sinee EE^ is 


strietly quasieoneave in Pk [54|; for given P^, the unique Pc ean be ealeulated from the left 


hand side of ( |E.15| ). 

The remaining task is to prove that Pc is strietly inereasing with P^. Define the left hand side 


of (E.151 as /(Pfc). Eunetion /(P^) is obviously differentiable in eaeh interval of the subregions, 
i.e., [gi, gi+]\ , / = 1, • • ■ , r^, whieh ean be ealeulated as 


fiPl) = -C"k{P;)Ck{P*k){Ck'{Pl))-^ > 0,Pfc e {gugi+i),l = I,-- - ,rfc. 


(E.16) 


where the inequality follows from the faets that Ck{Pk) is positive and C"k{Pk) is negative 


aeeording to (E.4) in Eemma 10. Combining with the faet that f{Pk) is a eontinuous funetion 
of P^, we eonelude that funetion /(P^) is a strietly inereasing funetion w.r.t. P^ and the proof 
for the first part of Theorem 4 eompletes. 

Einally, we prove the last part of Theorem 4: the EE deereases with the eireuit power 
eonsumption. Supposing P^^^ > P^\ define Pj:^'’ and Pjf'^ respeetively as the optimal solutions 
for given P^'^ and PqK Then, we have 


(iK_ 


P’ + Pc 


( 1 ) 


< 


H" + Pc 


( 2 ) 


< 


Ck{P, 


(2). 


= EEkiPg'^), 


PP^ + P 


( 2 ) 


(E.17) 


( 2 ) 

where the seeond inequality follows from the assumption that P^ is the optimal solution for 

( 2 ) 

given P^ . Erom the above inequalities, we know that the EE deereases with Pc. 
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